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Abstract 

Using numerical results obtained within two models describing vortex mat- 
ter (interacting elastic lines (Bose model) and uniformly frustrated XY-model) 
we establish universal properties of the melting transition within the linelike 
regime. These properties, which are captured correctly by both models, in- 
clude the scaling of the melting temperature with anisotropy and magnetic 
field, the effective line tension of vortices in the liquid regime, the latent heat, 
the entropy jump per entanglement length, and relative jump of Josephson 
energy at the transition as compared to the latent heat. The universal prop- 
erties can serve as experimental fingerprints of the linelike regime of melting. 
Comparison of the models allows us to establish boundaries of the linelike 
regime in temperature and magnetic field. 
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It is now well established that the vortex lattice in a clean superconductor melts in a first- 
order transition due to thermal fluctuations. The melting transition has been observed in 
a number of experiments using different techniques, including transportEHl and direct mea- 
surements of the magnetization jumpcrQ and the latent heatffl associated with the meltin E 
transition. 

Melting of the vortex lattice in real superconductors is a complicated phenomenon. Sev- 
eral different regimes exist depending on the parameters of the superconductors and the 
magnetic field. The simplest and most natural picture of the melting transition in three- 
dimensional superconductors is that a low temperature crystal of straight vortex lines trans- 
forms into a liquid of entangled lines.! In real layered superconductors this picture does not 
work either at too high or too low fields. At high fields, the vortex lattice is more adequatly 
described as a system of weakly interacting two-dimensional lattices. Melting in this situa- 
tion occurs in a pointlike way and is accompanied by evaporation of vortex lines. At low 
fields, on the other hand, the melting temperature approaches the fluctuation region and 
vortex degrees of freedom start to interact with other superconducting degrees of freedom 
(spin wave excitations and vortex loops). The contribution of these degrees of freedom to 
the latent heat of the melting starts to increase abruptly as the temperature aproaches T c 
possibly leading to a significant change in the nature of the melting transition.BH 

As a quantitative theory of vortex lattice melting has yet to be developed, a large amount 
of work has been invested into numerical simulations of the vortex system. A complete 
picture can be worked out by direct simulation of the Ginzburg-Landau functional, but 
this would require enormous computational effort. Rather a number of simplified models 
have been used, such as interacting elastic lines (or the Bose model)t3t3, the frustrated 
XY-mode]liac3, the Lattice London Model (or lines and ioops)EHI, and the Lowest Landau 
Level approximation.i3@ Being approximations, all these models have their limited regime 
of applicability where they can be expected to make reliable predictions for the melting 
transition and for the properties of different vortex phases. Unfortunately, it is difficult even 
to make reliable analytical predictions for the range of applicability of a specific model. In 
this work we therefore use a more pragmatic approach of comparing the numerical results 
for two rather different models: the Bose model and the XY-model. We find that the results 
obtained from these models agree over a wide range of fields and temperatures, where the 
melting transition can be described as linelike. We show that in the linelike regime the 
transition is characterized by a set of universal properties, which can serve as fingerprints of 
this regime. Comparison of the models also allows us to establish boundaries of the linelike 
regime. 



I. DESCRIPTION AND COMPARISON OF THE MODELS 

The phenomenology of superconductors is based on the Ginzburg-Landau model, which 
expresses the energy of superconductor via the magnetic induction, B, the modulus of 
the order parameter, and its phase, <fi = arg\l/. The model is completely defined 
by the London penetration depth A, the coherence length £, and the anisotropy factor 
7. Since the high-T c materials are strongly type II superconductors, i.e., k = A/£ ^> 
1, the modulus of the order parameter is suppressed only in the vicinity of vortex cores. 
This allows us to describe superconductors within the London approximation which neglects 
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fluctuations in the modulus of the order parameter everywhere except in the cores of the 
vortices. The approximation breaks down in the vicinity of the upper critical field H c2 and 
in the fluctuation region near the transition temperature T c . The London approximation is 
used by both models studied in this paper. The frustrated XY-model then makes two further 
approximations: It neglects the fluctuations in the magnetic induction, usually called the 
frozen-field or infinite-lambda approximation, and introduces a lattice on which the phase of 
the order parameter is defined. The model desribes layered superconductors only in terms 
of the phase distribution 0(n), defined on a three-dimensional grid n = (n x ,n y ,n z ). The 
energy functional of this model is given by (see, e.g., Ref. ESI). 



^[0(n)] = > ]{ J J2 y [<K n + d a ) - 0(n) - a Q (n)] - ~ 2 cos [0(n + d z ) - 0(n)] } , (1) 

a=x,y ^ 

The energy scale of the system is the phase stiffness, J, definied as 

= s$ 2 = seo 

7r(47rA a6 ) 2 7T ' 1 ' 

where s is the layer spacing and we have introduced the scale for the vortex line energy 
Eo = ($o/47rA a b) 2 with <3> being the flux quantum and \ a b is the penetration depth for 
supercurrents flowing in the afe-planes. The anisotropy is given by 7 = A c /A a fc, where 
A c is the penetration depth for supercurrents flowing the direction along the c-axis. The 
phase 0(n) is defined on a cubic lattice with the lattice spacing s, the dimensionless vector 
potential is defined as a = (0, 2nfn x /^ , 0), with / being the fraction of the lattice cells 
filled by vortices, and d Q are unit vectors. The phase interaction V{4>) is a 27r-periodic 
function with the Taylor expansion V(<f>) — V(0) ~ 2 /2 as — > 0. An obvious choice is to 
use V(<f>) = — cos(0), but this produces large barriers for the motion of the vortices between 
the lattice cells. A natural improvement is therefore to choose a function which minimizes 
these barriers, as has been done in Ref. |29| . 

The Bose model, on the other hand, describes the vortex system entirely in terms of 
the vortex degrees of freedom. A vortex is treated as an elastic string, interacting with 
the other lines through a screened Coulomb potential represented by the modified Bessel 
function K (-R/A a b). The free energy functional is then given byH 

JTB.W] = T <JV ? ( ^ V + T £ „K f £lH , (3) 




with L z is the thickness of the sample, Rj(z) is a two-dimensional vector describing the 
position of the vortex, and Si is elasticity of a vortex line. This single-line elasticity is an 
effective quantity, which we choose to reproduce the relevant tilt modulus of the vortex 
lattice. Two effects complicate the relation between this parameter and the parameters for 
superconductors. First, it is well known that the tilt energy of vortex line is nonlocal and 
the effective line stiffness for deformations with wave vector k z is proportinal to \n{^f/k z ^). 
The k z relevant for melting are given by k z ~ 7/aQ. This leads to the estimated 



7 2 \2y/n£ 
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where a is the lattice spacing of a triangular lattice, i.e., <&o/B = \^3al/2. The second 
effect becomes important in the temperature range close to the fluctuation region where 
phase and vortex fluctuations start to suppress the Josephson coupling. This suppression 
enhances the anisotropy and reduces the effective line tension of vortices E\ as compared to 
the estimate (^). In the following we introduce the effective anisotropy parameter e of the 
Bose model as e 2 = £i/eq, bearing in mind that e ~ I/7, with the main difference coming 
from the logarithm in Eq. (ffl). 

Both models have their limitations in describing real superconductors. The Bose model 
describes a continuous anisotropic system whereas the XY-model takes the layered structure 
of the superconductor into account. Thus one would not expect the Bose model to be 
applicable when the typical wavelength in the z-direction becomes comparable to the layer 
spacing. As will be shown below, the correlation length in the ^-direction for the Bose model 
is l z = aao/7, where a ~ 6 is a numerical constant!^ Therefore, we expect the Bose model 
to be invalid when ao < 7s/a or for B <La 2 B cr . Numerical results show that line melting 
works at least for fields as large as 10-B cr .c3 For larger fields, the Josephson coupling between 
the layers is strongly suppressed by thermal fluctuations, invalidating the use of a simple 
line model. 

Both models neglect the electromagnetic coupling between the superconducting layers. 
This coupling is relevant for very anisotropic materials when 7s > ao, A.0 The XY-model 
further neglects the screening of the vortex interaction due to the magnetic field, which 
becomes relevant close to T c when ao ^ A (T) , and therefore does not describe the reentrant 
behavior of the melting line which is expected for low fields. 

When the melting temperature approaches the fluctuation region thermally activated 
vortex loops start to influence the thermodynamic properties of superconductor strongly. 
The Bose model only contains the field induced vortices as degrees of freedom and neglects 
this effect. The XY-model allows for both Gaussian phase fluctuations and thermally in- 
duced vortex loops in addition to the field induced vortices. The extra degrees of freedom 
give a large contribution to the specific heat of the system: Typical numerical data from the 
XY-model show a specific heat which grows with temperature and has a broad maximum 
at T c . On too of this specific heat there is a small sharp peak due to the first-order melting 
transition.0'cl These results agree well with experimental observations.0! The important is- 
sue is to understand how the field-induced vortices couple to the other fluctuations. Outside 
the fluctuation region, the coupling only leads to a weak renormalization of vortex interac- 
tion and the effective line energy, as is apparent from a number of observations: To begin 
with, the shape of the melting line is well reproduced by the Lindemann criterion, which 
only includes field induced vortices with unrenormalized interactions (see, e.g., Refs. |32] and 
0). Furthermore, it has recently been shown that the field and temperature dependence of 
the discontinuity at the vortex lattice melting transition also can be estimated using a sim- 
ple line model and the mean-field temperature dependence of Eo(T)M. On the other hand, 
recent simulations have shown that if melting takes place in the vicinity of the fluctuation 
region, then a liquid phase immediatly above the transition is not simply a liquid of inter- 
acting elastic lines, as the Bose model assumes, but rather an infinite interconnected cluster, 
in which individual lines loose their identities.!^ Therefore it is natural to test the adequacy 
of the line picture and establish boundaries of its applicability by direct comparison of the 
numerical results of the Bose model with the more microscopic XY-model. 
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II. UNIVERSAL PROPERTIES. COMPARISON OF NUMERICAL RESULTS 



Important conclusions can be drawn already from the relevant scales for energy, length, 
and magnetic field in the expressions ([I]) and (13]). For the XY- model, the relevant scale for 
the magnetic field is the crossover field B CT = $ /(7 s ) 2 an d the relevant energy scale is the 
phase stiffness jM Thus, melting lines for systems with different anisotropy factors should 
collapse on to one line as long as the field is measured in units of B CT and temperature is 
measured in units of J, as has indeed been demonstrated experimentally.E3ia In the region of 
small fields B < B CT and beyond the fluctuation region T < J the lattice is expected to melt 
in a linelike fashion, i.e., the lines are expected to retain their identity above the transition 
and the Josephson coupling between the layers is only weakly suppressed. In this regime 
one can expand the Josephson coupling energy with respect to interlayer phase difference 
and make a continuous approximation, which leads to an even simpler scaling property: The 
melting temperature should scale as 

T m = A m Jy/B cr /B. (5) 

Simulations show that this scaling indeed works very well for fields B cr ;$ B ;$ 105 cr with 
A-rn ~ 0.33 even though a noticable suppression of the Josephson coupling is observed near 
T m (down to 64% of the bare coupling) .H The scaling form (H) assumes, in fact, that only 
vortex degrees of freedom participate in the melting transition and regular phase fluctuations 
play negligible role at T m . This approximation breaks down when the melting temperature 
approaches the fluctuation region, which corresponds to T m > J or B < B cr . Regular 
phase fluctuations and thermally activated vortex loops suppress both the in-plane phase 
stiffness and the Josephson interlayer coupling. This leads to a renormalization of vortex 
interactions and their tilt stiffness. Quantitatively these effects can be characterized by 
the helicity moduli T x and T 2 in the Meissner state. A natural generalization of the scaling 
form (|5|), which takes into account this renormalization, can be obtained by the replacements 
J — > T x and 7 — > with 7y = T z /T x being the anisotropy of the helicity moduli. This 
leads to a generalized scaling relations for melting temperature 



T m = A m T x ^%/B{ lY s) 2 (6) 

For most of the phase diagram melting occurs in the region where regular phase fluctuations 
are weak and one can use formulas for T Xj2 with small fluctuation corrections (see Appendix) 



| ^ yWT / ln(32 7 2 ) - 1 
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4J V* 27r 7 2 
T(ln(32 7 2 ) - 1) 
4^J 



(7) 
(8) 



with V^' = —d A V{0)/d0 4: at 9 = and V(9) is the phase interaction function from Eq. (|T]). 
To check the scaling relation (^|) we used all availible data on the melting transition for XY 
models with different anisotropies and filling factors .SHEOH Fig. [1] shows the dependence 
of T m /T x vs l/y/fjY, which according to Eq. (||) should be a straight line. We found 
that in a suprisingly wide range of fields and anisotropies, including a substantial part of 
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the fluctuation region, the scaling works very well with A m « 0.4. At low fields the 
generalized scaling relation extends down to B « 0.1-B cr . 

Near the transition point T c = a c J the helicity moduli vanish according to the XY-scaling 
laws 



rp\ 2v 



T.« J [«„--) , (9) 



2u 



with v m 1/3. Note that above representations include the possibility of a temperature 
dependent J near the mean field transition temperature T c0 , J = J (T c0 — T)/T c0 . Assuming 
that the scaling law (H) can be extended to the fluctuation region we conclude that in the 
vicinity of the transition temperature T c the melting field should scale as B m oc (T c — T) 4//3 . 
This scaling indeed describes the behavior of experimental melting line better than the 
"mean field" scaling B m oc (T c q — T) 2 , which follows from scaling relation ([5]) (see, Refs. 



The Bose model is able to describe the melting transition for arbitrary ratio clq/\. There 
is no general scaling relation which holds for arbitrary values of this ratio. However, in the 
region of the high vortex density a <C A, where the interaction can be approximated with a 
logarithm, the free energy only depends on the dimensionless parameter A = T/ao\/2£i£o-0 
In particular, the melting takes place at A m = 0.06220 giving the following scaling relation 
for the melting temperature. 

T m « 0.088ee a . (11) 

This scaling form corresponds to the low field scaling of the XY-model given by Eqs. fl5|) 
and @. Rewriting Eq. (|5]) in terms of a we obtain 



m 0.33 / VOEnan 

T m « \/ -V— « 0.098e a /7- (12) 

7T y 2 7 

Thus, the two models do not only agree on the shape of the melting line, but they also agree 
quantitatively on the position of the melting transition. The result is consistent with the 
assumption that, according to Eq. (Q), we have e > I/7 and we find 

e« I.H/7. (13) 

For finite range interactions the melting temperature is lowered; simulations with A = 
1.06or) give T m rs 0.084e£:oOo, i-e., the melting temperature is lowered by approximately 
5%B 

We now compare properties of the line liquid near the transition. Both experiments and 
simulations show the vortex liquid to be heavily entangled so that the phase coherence in 
the direction parallel to the applied field is destroyed by melting. The entanglement can be 
quantified in terms of the transverse wandering of the vortex lines, 

w(z) = ([R(z)-R(0)} 2 )/al (14) 
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It is easy to show that we have 



2T z 

w(z) = —-, (15) 
Eido ao 

for non-interacting lines. Since interactions between the lines are unimportant for short 
distances, this is also the result for small z in the Bose model. For larger z, the wandering is 
suppressed due to the interaction between the vortices but we still expect w(z) oc zj. The 
numerical result is 

2T z 

w{z) = 0.27 , (16) 

£;Oo do 

implying that the line stiffness is increased by a factor 3.7 due to the interaction with 
other lines in the liquid. For the XY-model as well as for real superconductors the linear 



dependence ([15]) is only approximate due to nonlocality of the line tension. In Fig. |2| we 
show w(z) for two different values of 7 in the XY-model, just above the melting transition. 
We also show w(z) from the Bose model, where we have chosen e to give the best fit. This 
yields, 

w(z) » lAQz/ea 1.46/7, (17) 

which again is in reasonable agreement with with Eq. (Pi). For comparison the plot also 
shows w(z) in the lattice just below the melting transition. One can see that line wandering 
agrees very well for the two models also in the lattice state. In particular an, initial growth 
of w(z) is similar to the liquid state and it starts to saturate at approximately the same 
length scale along z-direction. For the Bose model no entanglement is observed in crystal 
state and line wandering saturates at w(z) m 0.15. For the XY-model w(z) continues to 
grow due to rare reconnection events, with slope much smaller than in the liquid state. This 
small slope increases with decreasing anisotropy. 

It is convenient to have a single number characterizing the line wandering in the vortex 
liquid. A suitable quantity is the entanglement length l e , which is defined as a length at 
which average line displacement becomes equal to the radius of the Bravais cell, 

w(/e) = 7^ w °- 276 - ( 18 ) 

From Fig. |2] we find an entanglement length l e « 6.600/7, implying that the vortex system 
can become fully entangled already in a thin sample. 

We now turn to the thermodynamic properties of the melting transition obtained from 
simulations. A main characteristic of the melting transition is the latent heat or entropy 
jump. As is now generally accepted, it is important to take the internal temperature de- 
pendence of the free energy functional T into account when computing the energy of a 
superconductor using Ginzburg-Landau theoryBii As a result, the energy is not just the 
average of T but is rather given by 

E =^)- T {w)> (19) 
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where (...) indicates the Monte Carlo average. In the XY- model, the temperature depen- 
dence of T is only due to the temperature of the energy scale J = ssq/ti. The same is true 
for the Bose model if we assume a logarithmic interaction between the vortices, i.e., when 
A > a . We then obtain the simple result 

E=(l (20) 



e dT 

This equation, together with numerical data for the discontinuity in (J 7 ), produces results 
which are consistent with experiments.il 

It should be remembered, however, that the statistical weight of a particular vortex 
configuration is still given by T. It therefore makes sense to study the configurational energy 
E c = (J 7 ) and the corresponding entropy, S c = E c /T = (JF)/T. It is this configurational 
entropy which we expect to be proportional to the number of degrees of freedom of the 
vortex system. In order to be able to compare the results from the two models, we only 
consider the energy per vortex and length, 

L ' 1 (21) 



NL 7 NL, 



The Bose model giveail 



Ae c « 0.013e (22) 

and the XY-model typically produces a jumpH 

Ae c ^ 0.016e - (23) 

It is interesting to note that the jump in the configurational entropy per entanglement length 
at the transition is 

As c l e = Ae c l e /T m w (0.9 - l.l)k B . (24) 

Thus, loosely speaking, the system has gained one degree of freedom per field-induced vortex 
and entanglement length. The result is again perfectly consistent with linelike melting. 

The discontinuity in magnetization can be computed from the entropy jump using the 
Clausius-Clapeyron relation, 

4vr dT ' v ; 

where H m (T) is the applied external field at melting. If we use the approximation H m {T) w 
B(T), which again is true for A > ao, we can rewrite this as 

^ = 8ttA 2 A P , (26) 

where p = B/<&q is the density of the flux lines. If we further assume the latent heat to be 
independent of the interaction range, we obtain a simple expression for the magnetization 
jump, 
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AB » 5.2 x KT 4 $ /A. (27) 

This result explains why the magnetization jump vanishes in the limit A — > oo and is in 
good agreement with experimental results for YBCO. An advantage with the Bose model 
is that we can do a simulation with finite A, allowing the density to fluctuate, and test the 
consistency of the simulation with Clausius-Clapeyron directly.EJ 

It is also interesting to examine what fraction of the latent heat comes from the jump 
in tilting energy of vortices which corresponds to the Josephson coupling between the layers 
in the XY model. It was shown that for the Bose model the jump in the aver ag e tilt 
energy Aeuu is exactly half of the latent heat Ae c , where euu = Ei((dH/ dz) 2 ) /2.li3 This 
can also be considered as an indication of the linelike regime. To investigate this property 
for the XY model one has to study the relative jump of the Josephson energy Aej/Ae c 
with ej = — / J/(7 2 l(cos[0(n + d z ) — <p(n)]). Simulations shows that for a wide range of 
fields Aej ~ 0.5Ae c xB Moreover, using the fact that the free energy must be continuous at 
melting, the following relation can be derived 

Aej = <9mT m 

Ae c <91n 7 2 ' 1 ' 

which shows that this property is a consequence of the scaling relation fl5|). 

To conclude, we have compared and discussed numerical result from two widely used 
models for the vortex system, the XY-model and the Bose model. We have shown that 
the models agree even quantitatively over a large part of the phase diagram, where we 
conclude the melting to be linelike. We have identified a number of universal properties, the 
entanglement length and configurational entropy per entanglement length, which we believe 
are generic to the transition in the linelike regime. 
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The HPCRF is funded principally by the U.S. Department of Energy Office of Scientific 
Computing. 



III. APPENDIX: FLUCTUATION CORRECTIONS TO HELICITY MODULI IN 

XY-MODEL 

The helicity moduli T x and T 2 describe the supercurrent responses to phase rotations. 
In the Meissner state they are determined by the following thermodynamic averages 

J 2 l\ 

T X = J (V" (0(n + d x ) - 0(n))> - — - ( £V (0(n + d.) - 0(n)) 



T 2 = 4(cos(0(n + d,)-0(n))> 




(29) 

^sin(0(n + d 2 )-0(n)) ), (30) 




9 



where N is the total number of grid sites, V = dV (<fi) j 'dcj) , V" = d 2 V {<p) / d<p 2 , and (. . .) 
notates average with the energy functional (1). At finite temperatures the helicity moduli 
are suppressed by fluctuations below their zero temperature values J and J/7 2 



T, 



J(l 

J_ 
-2 



1 , 



<(0(n + d x .)-0(n)) 2 > 
(0(n + d z )-0(n)) 2 > 



(31) 
(32) 



where = d A V (<p) / dcf)^^ . Phase fluctuations at small temperatures are determined by 
the Gaussian energy functional which is obtained by expansion of Eq.(l) 



J 



n \_a=x,y 

Using this functional we calculate 



(33) 



n + 4)-0(n)) 2 > 
<(0(n + 4)-0(n)) 2 > 



T 
2J 
T 
2^J 



1 - 



27T7 2 

(ln32 7 2 - l) 



(ln32 7 2 - l) 



(34) 
(35) 



Substituting these expressions into Eqs. (|3"ID and (|32| ) we obtain Eqs. (7) and (8) of the 
paper. 
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FIG. 1. Scaling plot of melting temperature in the line melting regime for XY models with 
different anisotropy factors. Data are taken from Refs. 2C, 2^,21, 23. Linear fits are made for data 



from Refs. 29.23 
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FIG. 2. Comparison of line wandering in the crystal and liquid states near the melting point 
for the Bose model and the frustrated XY-model. The two models show very similar behavior both 
with regards to the asymptotic and the initial slope in the liquid. Likewise, the saturation of the 
line wandering in the crystal is captured correctly by both models. 
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